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Quasi-equilibrium models for aggregate variables are widely-used throughout finance and eco- 
nomics. The validity of such models depends crucially upon assuming that the systems' participants 
behave both independently and in a Markovian fashion. 

We present a simplified market model to demonstrate that herding effects between agents can 
cause a transition to boom-and-bust dynamics at realistic parameter values. The model can also be 
viewed as a novel stochastic particle system with switching and reinjection. 

PACS numbers: 89.65. Gh, 89.75.Fb, 89.75.Da 



I. INTRODUCTION 

In the physical sciences using a stochastic differential 
equation (SDE) to model the effect of exogenous noise 
upon an underlying ODE is often straightforward. The 
noise consists of many uncorrelated effects whose cumu- 
lative effect is well-approximated by a Brownian process 
B s , s > and the ODE df = a(f, t) dt is replaced by an 
SDE df = a(f, t) dt + b(f, t) dB t . 

However, in financial and socio-economic systems the 
inclusion of exogenous noise (ie new information enter- 
ing the system) is more problematic — even if the noise 
can be legitimately modeled as a Brownian process. This 
is because such systems are themselves the aggregation 
of many individuals or trading entities (referred to as 
agents) who typically 

a) have differing interpretations of the new information, 

b) act differently depending upon their own recent his- 
tory (ie non-Markovian behaviour), and 

c) may not act independently of each other. 

The standard approach in neoclassical economics and 
modern finance is simply to 'average away' these awk- 
ward effects by assuming the existence of a single repre- 
sentative agent as in macroeconomics pQ, or by assum- 
ing that the averaged reaction to new information is cor- 
rect/rational, as in microeconomics and finance [SJ [3J. In 
both cases, the possibility of significant endogenous dy- 
namics is removed from the models resulting in unique, 
Markovian, (quasi)-equilibrium solutions. 

In reaction to this, many Heterogeneous Agent Models 
(HAMs) have been developed |1] that simulate the agents 
directly. These have demonstrated that it is relatively 
easy to generate aggregate output data, such as the price 
of a traded asset, that approximate reality better than 
the standard averaging-type models. In particular the 
seemingly universal 'stylized facts' of financial mar- 
kets such as heteroskedasticity (volatility clustering) and 
leptokurtosis (fat-tailed price-return distributions result- 
ing from booms-and-busts) have been frequently repro- 
duced. However, the effects of such research upon main- 
stream modeling have been minimal perhaps, in part, 
because some HAMs require fine tuning of important pa- 
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FIG. 1. All the M signed particles are subject to a horizontal 
stochastic forcing and they also diffuse independently. Minor- 
ity particles also drift downwards at a rate proportional to the 
imbalance. When a particle hits the boundary it is reinjected 
with the opposite sign and a kick is added to the bulk forcing 
that can trigger a cascade (see text). 



rameters, others are too complicated to analyze, and the 
plethora of different HAMs means that many are mutu- 
ally incompatible. 

The purpose of this report is to investigate the range 
of validity of the quasi-equilibrium solutions obtained by 
ignoring endogenous effects such as a), b) and c) above. 
We do this by introducing a simplified version of the mod- 
eling framework introduced in [3 [8] that can also be de- 
scribed as a particle system in two dimensions (Figure!!]). 



II. A STOCHASTIC PARTICLE SYSTEM WITH 
REINJECTION AND SWITCHING 

We define the open set D C K 2 by D — {(x, y) : —y < 
x < y, y > 0. There are M signed particles (with states 
+1 or —1) that move within D subject to three different 
motions. Firstly there is a bulk Brownian forcing B t in 
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the x-direction that acts upon every particle. Secondly, 
each particle has its own independent two-dimensional 
diffusion process. Thirdly, for agents in the minority 
state only, there is a downward (negative y-direction) 
drift that is proportional to the imbalance. 

When a particle hits the boundary dD it is reinjected 
into D with the opposite sign according to some prede- 
fined probability measure. Finally, when a particle does 
switch the position of the other particles is kicked in the 
x-direction by a (small) amount ±§y, t > 0, where the 
kick is positive if the switching particle goes from the —1 
state to +1 and negative if the switch is in the opposite 
direction. Note that the particles do not interact locally 
or collide with one another [5]. 



A. Financial market interpretation 

We take as our starting point the standard geometric 
Brownian motion (gBm) model of an asset price p t at 
time t with po = 1. It is more convenient to use the log- 
price r t = lnp t which for constant drift a and volatility 
b is given by the solution r t = at + bB t to the SDE 



state = -1 



dr t = a dt + b dB t 



(1) 



Note that the solution r t depends only upon the value 
of the exogenous Brownian process B t at time t and not 
upon {-E> s }s = o- This seemingly trivial observation im- 
plies that r t is Markovian and consistent with various 
notions of market efficiency. Thus gBm can be consid- 
ered a paradigm for all economic and financial models in 
which the aggregate variables are in a quasi-equilibrium 
reacting adiabatically to new information. 

The instantaneous translation of new information B t 
into price changes is effected by 'fast' agents who will not 
be modeled directly. However, we posit the existence of 
M 'slow' agents who are primarily motivated by price 
changes rather than new information and act over much 
longer timescales (weeks or months). At time t the i th 
slow agent is either in state Si(t) = +1 (owning the asset) 
or Si(t) = — 1 (not owning the asset) and the sentiment 

a(t) e [-1, 1] is defined as a(t) = J2Zi s *(*)- The i& slow 
agent is deemed to have an evolving strategy that at time 
t consists of an open interval (Li(t), Ui(t)) containing the 
current log-price r* (see Figure^. The i th agent switches 
state whenever the price crosses either threshold, ie rt = 
Li{t) or Ui(t), and a new strategy interval is generated 
straddling the current price. 

We assume in addition that each threshold for every 
slow agent has its own independent diffusion with rate 
oti (corresponding to slow agents' independently evolving 
strategies) and those in the minority (whose state differs 
from sgn(er)) also have their lower and upper thresholds 
drift inwards each at a rate Ci\a\, Ci > 0. 

These herding constants Ci are crucial as they provide 
the only (global) coupling between agents. The inward 
drift of the minority agents' strategies makes them more 
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FIG. 2. A representation of the model showing two agents in 
opposite states. Agent i is in the +1 state and is represented 
by the two circles at Li(t) and Ui(t) while agent j is in the 
— 1 state and is represented by the two crosses. 



likely to switch to join the majority. Herding, and other 
mimetic effects, appear to be a common feature of finan- 
cial and economic systems. Some causes are irrationally 
human while others may be rational responses by, for ex- 
ample, fund managers not wishing to deviate too far from 
the majority opinion and thereby risk severely underper- 
forming their average peer performance. The reader is 
directed to [5] for a more detailed discussion of these and 
other modeling issues. 

Finally, changes in the sentiment a feed back into the 
asset price so that gBm ([I]) is replaced with 



dr t = a dt + b dB t + kAct 



(2) 



where k > and the ratio n/b is a measure of the relative 
impact upon r t of exogenous information versus endoge- 
nous dynamics. Without loss of generality we let a = 
and b — 1 by setting the risk-free interest rate to zero 
and rescaling time. 

One does not need to assume that all the slow agents 
are of equal size, have equal strategy-diffusion, and equal 
herding propensities. But if one does set oti = a and 
Ci — C Vi then one obtains the particle system simply 
by defining the position of the i th particle as (xi,yi) = 
( Ui + Li — r t , Ui ~ Li ), In other words, the bulk stochastic 
motion is due to exogenous noise; the individual diffu- 
sions are caused by strategy-shifting of the slow agents; 
the downward drift of minority agents is due to herd- 
ing effects; the reinjection and switching are the agents 
changing investment position; and the kicks that occur at 
switches are due to the change in sentiment affecting the 
asset price via a linear supply /demand-price assumption. 



B. Limiting values of the parameters 

There arc different parameter limits that are poten- 
tially of interest. 

1) M —¥ oo In the continuum limit the particles are re- 
placed by a pair of evolving density functions p + (x, y, t) 
and p~{x 1 y,t) representing the density of each agent 
state on D — such a mesoscopic Fokker-Planck descrip- 
tion of a related, but simpler, market model can be found 
in |10) . The presence of nonstandard boundary condi- 
tions, global coupling, and bulk stochastic motion present 
formidable analytic challenges for even the most basic 
questions of existence and uniqueness of solutions. How- 
ever, numerical simulations strongly suggest that, minor 
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discretization effects aside, the behaviour of the system 
is independent of M for M > 1000. 

2) B t — > As the external information stream is reduced 
the system settles into a state where a is close to cither 
±1. Therefore this potentially useful simplification is not 
available to us. 

3) a — > or oo In the limit a — ¥ the particles do not dif- 
fuse ie. the agents do not alter their thresholds between 
trades/switches. This case was examined in [IT] and the 
lack of diffusion does not significantly change the boom- 
bust behaviour shown below. On the other hand, for 
a ^> max(l,C) the diffusion dominates both the exoge- 
nous forcing and the herding/drifting and equilibrium- 
type dynamics is re-established. This case is unlikely in 
practice since slow agents will alter their strategies more 
slowly than changes in the price of the asset. 

4) C — » This limit is the focus of the report. When 
C = the particles are uncoupled and if the system 
is started with approximately equal distributions of ±1 
states then a remains close to 0. Thus ^ reduces to 
and the particle system becomes a standard equilib- 
rium model — agents have differing expectations about 
the future which causes them to trade but on average the 
price remains 'efficient' []. In Section III we shall observe 



that endogenous dynamics arise as C is increased and the 
equilibrium solution is no longer stable. 
5) k — y For k > one agent switching can cause an 
avalanche of similar switches, especially when the sys- 
tem is highly one-sided with |<r| close to 1. When re = 
the particles no longer provides kicks (or affect the price) 
when they switch although they are still coupled via 
C > 0. The sentiment a can still drift between —1 and +1 
over long timescales but switching avalanches and large, 
sudden, price changes do not occur. 



III. PARAMETER ESTIMATION, NUMERICAL 
SIMULATIONS AND INSTABILITY 



In all the simulations below we use M = 1000 and 
discretize using a timestep h — 0.000004 which corre- 
sponds to approximately 1/10 of a trading day if one 
assumes a daily standard deviation in prices of « 0.6% 
due to new information. The price changes of 10 con- 
secutive timesteps are then summed to give daily price 
return data making the difference between synchronous 
vs asynchronous updating relatively unimportant. 

We choose a = 0.2 so that slow agents' strategies dif- 
fuse less strongly than the price does. A conservative 
choice of re = 0.2 means that the difference in price be- 
tween neutral {a — 0) and polarized markets a = ±1 is, 
from exp(0.2) 22%. 

After switching, an agent's thresholds are chosen ran- 
domly from a Uniform distribution to be within 5% and 
25% higher and lower than the current price. This allows 
us to estimate k by supposing that in a moderately polar- 
ized market with |er| = 0.5 a typical minority agent (out- 
numbered 3-1) would switch due to herding pressure af- 




FIG. 3. In each picture, red corresponds to output with the 
parameters used in the text while blue represents outputs for 
the gBm pricing model (|TJ with the same B B . The top left 
figure shows the prices pit); top right is sentiment vs time; 
bottom left plots the daily price changes; and bottom right 
shows the cumulative log-log plot of daily price changes that 
exceed a given percentage. 



ter approximately 80 trading days (or 3 months, a typical 
reporting period for investment performance) [12 . The 
calculation 80C|er| = | ln(0.85)|/0.00004 gives C « 100. 
Finally, we note that no fine-tuning of the parameters is 
required for the observations below. 

Figure [3] shows the results of a typical simulation, 
started close to equilibrium with agents' states equally 
mixed and run for 40 years. The difference in price his- 
tory between the above parameters and the equilibrium 
gBm solution is shown in the top left. The sudden mar- 
ket reversals and over-reactions can be seen more clearly 
in the top right plot where the market sentiment under- 
goes sudden shifts due to switching cascades. These re- 
sult in price returns (bottom left) that could quite easily 
bankrupt anyone using excessive financial leverage and 
gBm as an asset pricing model! Finally in the bottom 
right the number of days on which the magnitude of the 
price change exceeds a given percentage is plotted on log- 
log axes. It should be emphasized that this is a simplified 
version of the market model in [5] and an extra parame- 
ter that improves the statistical agreement with real price 
data (by inducing volatility clustering) has been ignored. 

To conclude we examine the stability of the equilib- 
rium gBm solution using the herding parameter C as 
a bifurcation parameter. In order to quantify the level 
of disequilibrium in the system we record the maximum 
value of \a\ ignoring the first 10 years of the simulation 
(to remove any possible transient effects caused by the 
initial conditions) and average over 20 runs each for val- 
ues of < C < 40. All the other parameters and the 
initial conditions are kept unchanged. 
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The results in Figure [4] show that for values of C as 
low as 20 the deviations from the equilibrium solution 
are as large as the system will allow, with the large ma- 
jority of agents being in the same state at some point 
during the simulation. This is lower than the value of 
C = 100 estimated above by a significant margin and it 
should be noted that there are other phenomena, such 
as new investors and money entering the asset market 
after a bubble has started, and localized interactions be- 
tween certain subsets of agents that might cause further 
dcstabilization. 
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FIG. 4. A measure of disequilibrium |(j| max averaged over 20 
runs as the herding parameter C changes. 



IV. CONCLUSIONS 

Financial and economic systems are subject to many 
different kinds of inter-dependence between agents and 



potential positive feedbacks. However, even those main- 
stream models that attempt to quantify such effects [T2] 
assume that the result will be a shift of the equilibria to 
nearby values without qualitatively changing the nature 
of the system. However we have demonstrated that at 
least one such form of coupling (herding) results in dis- 
equilibrium. Furthermore the new dynamics occurs at 
realistic parameters and is clearly recognizable as 'boom- 
and-bust'. It is characterized by long periods of low-level 
endogenous activity (long enough, certainly, to convince 
equilibrium-believers that the system is behaving adia- 
batically) followed by large, sudden, reversals involving 
cascades of switching agents triggered by price changes. 

The model presented here is compatible with existing 
(non-mathematized) critiques of equilibrium theory by 
Minsky and Soros [T51 1 14). Furthermore, work on related 
models to appear elsewhere shows that positive feedbacks 
can result in similar non-equilibrium dynamics in more 
general micro and macro-economic situations. 

Finally, the model has interesting links to other areas 
of mathematics and physics. In [8] it was shown that 
the switching cascades can be described using Queueing 
theory, with the price changes being equivalent to the 
busy-period of a queue. The model also shares similari- 
ties with other self-organizing systems such as the OFC 
earthquake model[TH]. And if one considers agents to be 
evolving hysteretic operators then results concerning the 
interaction of stochastic processes and hysteresis [16] and 
phase transitions may provide valuable insights into the 
price dynamics. 

The author thanks Michael Grinfeld and Rod Cross for 
numerous enlightening conversations and Julian Todd for 
writing the browser-accessible simulation of the particle 
system. The author also thanks Arjun Sanghvi for his 
additional computations and the NSF for his support. 
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